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COMPARISON OF THE »TBD” TOKEN 
MODEL TO AVAILABLE 
COMPETITORS 


CFD code to model aerosol flows on a large scale. Whereas DNS is the proper way to calculate 
flows on the scale of turbulent eddies, computers are still a long way away from being powerful 
enough to calculate flows in practical industrial and environmental systems. 

We have been developing a tool for the calculation of aerosol flows with the focus on the speed 
to reach the steady-state, efficiency and numerical accuracy. 
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KEY ELEMENTS OF "TBD" TOKEN MODEL 

• Basic Fluid Flow (1-D, 2-D, 3-D, Viscous, Compressible) 

• Internal \ External (e.g. pipes, ducts and rooms \ airplanes, cars) 

• Incompressible \ Compressible \ Transonic \ Supersonic \ Hypersonic (Mach 0-5) 

• Laminar \ Turbulent (e.g. small to large volumes) 

• Strong Coupling to Aerosol Dynamics & Chemistry (e.g. multi -physics) 

• Multi-Phase Aerosol Behavior (1-D, 2-D, 3-D, Polydisperse) 

• Aerosol Transport 

• Aerosol Dynamics 

• Chemically Reaction (1-D, 2-D, 3-D, Multi Species) 

• Multi Species Transport (Laminar \ Turbulent) 

• Strong Coupling to Aerosol Dynamics & Chemistry 

• Fluid phase and Surface Reactions 


The basic elements of the model are shown here. In order to simulate these a wide variety of 
cases, the coupled model has been made as general as possible. The model consists three 
modules ... 
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PROVEN TBD TOKEN APPLICATIONS 

• Prediction/Minimazation of Atmospheric Deposition 

• Prediction/Minimazation of Heat Exchanger Fouling 

• Aerosol Control in Nuclear Reactors 

• Design of Particle Characterization Equipment 

• Design of Dry Particle Inhalers 

• Design of Nebulizer Pulmonary Drug Delivery Systems 

• Design of Photocatalytic and Optical Particle Synthesis Processes 

• Design of High Surface Area Metal Catalyst Particle Synthesis Processes Design of 
Nanostructured Doped Silica Particles for Fiber Optic Preform Manufacture. 

• Design of Gas Phase Catalyst Synthesis Processes 

• Design of Carbon Nanotubes and Nano Onions Synthesis Processes 

• Design of Nanostructured Inhalable Drug and Drug Carrier Particle Synthesis 
Processes 

• Design of Nano-drug Particle Synthesis Processes 

• Design of Nano-particle Deposition Process 

• Etc. 


We have chosen to address two very different kinds of problems with this model to exemplify its 
generality. A brief background of these cases will be discussed... Aerosol Deposition in 
Combined Cycle Gas Turbines, and Aerosol formation and evolution in HSCT Nozzles and 
Plumes. 
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For example, the governing equations for a fluid flow can be written as a function of implicit 
primitive variables p, p, H and V as well as explicit variables p, y, k and MW. Chemical reactions 
and aerosol behavior can alter these quantities through various mechanisms. 

With regard to gas phase, for instance, reactions can absorb or release thermal energy thus 
changing enthalpy and pressure or altered species concentrations can change the p , MW, and 
p. A change in any such parameter can thus effect on the overall flow. 

Additionally, the presence of aerosols in a flow can create strong coupling effects. The 
formation of aerosols through condensation and nucleation of species, for instance, can release 
thermal energy, effecting state variables, OR large number of volume concentrations of aerosols 
can effect the effective viscosity of the fluid. 

Coupling in the reverse direction can also occur. Expansion of flow can alter saturation ratios 
leading to formation of aerosols, changing species concentrations and drive chemical reactions. 
Obviously these coupling effects can become quite involved. In two and three dimensions, 
additional factors due to spatial variation of parameters complicate the process even further. 
This presentation details and applies a method of modeling this intricate coupled behavior. 
Solutions presented today will focus on the coupling of fluid and gas phase species 
concentrations on aerosol behavior where reverse coupling is weak. When complete, the 
model will be able to model coupling in detail. 

Tow problems where this kind of coupling is important, namely: 

Aerosol deposition in Gas Turbines, and Pollutant Formation in FISCT Nozzles and Plumes 
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Creating a solution to the governing equations as analogous to getting water from a spring in 
the mountains to a pool in the a valley below. The water's position at the begging is like our 
initial conditions (out initial guess for the solution) and the position at the pool is our converged 
condition (where all the equations are satisfied). High in the mountains, our residual (error) is 
high. At the pool our residual is zero. 

The first thing to note is that there are many possible ways for the water to get from the spring 
to the pool. It can flow in a river or be carried in a helicopter. These represent alternate 
solution paths. There are also many possible springs that can reach the pool. These represent 
alternate initial conditions. 

Also note that all we care about is getting to the pool. The trip through the forest is not the 
important thing. If we could carry the water in a helicopter from the spring to the pool, that is 
fine. The problem is, we have to know where it is. 
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A coupled 2-D / Axisymmetric / 3-D flow and moment model incorporating the effects of 
nucleation, condensation, coagulation, diffusion, inertial impaction and thermophoresis has 
been developed. Computations have been presented with and RNS / aerosol model in general 
non-orthogonal coordinates for a variety of flows. The method has been shown to provide 
meaningful data on the complex behavior of aerosols under these flow conditions. 

A sample calculation of NASA's Variable Diameter Inlet geometry was used to 
demonstrate the capabilities of the flow module. In this complex, unsteady, turbulent 
geometry, the RNS methodology is found to accurately compute separation regions, normal and 
oblique shocks for internal \ external flow regions. The capabilities of the aerosol transport 
module have been demonstrated in a 3-D 180° turnaround duct. Complex flow-aerosol 
interactions are detailed. An analysis of gas to particle conversion in an incompressible to 
supersonic converging diverging nozzle demonstrates the detailed particle dynamics that can be 
modeled using this method. Preliminary results indicate the importance of multidimensionality 
in determining aerosol behavior in such geometries. Discrepancies between average or 
integrated values of aerosol properties in 2-D calculations can be as much as 10% as compared 
to 1-D calculations. Initial results indicate that the effects of aerosol dynamics on flow 
conditions are small for the relatively low conversion rates investigated. Validation of these 
results is ongoing. Analysis aerosol behavior shows the importance of correctly modeling all 
aspects of gas to particle conversion. For instance, it is shown that neglecting the nucleation 
process, even when seed particles are present, can result is large under-prediction of resultant 
aerosol volume concentrations and erroneous size distributions. A detailed analysis of aerosol 
effects on the flow is ongoing. 



It has been demonstrated that the present moment model is effective and efficient in 
solving for the complex behavior of polydisperse aerosols in general flow applications. Many 
important aerosol properties, such as size distributions, number and volume concentrations and 
depositional properties, can be calculated directly without having to sum up results for 
monodisperse aerosols over the entire particle spectrum. Research is presently underway to 
perform coupled flow and aerosol simulations under turbulent conditions in nozzles, plumes 
and inlets. 

Modelling 
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SCHEMATIC OF GRID SCALE OF DISPERSED PHASE 
PROCESSES OF INTEREST 


VOF, Embedded Interface 
& Level Set Methods 


Particle Tracking (Lagrangian) 
Sectional & Moment Methods 


Multi-Fluid 


Dynamic mechanisms such as coagulation, nucleation and condensation can become important 
in many instances. These mechanisms help define the chemical composition of the aerosols as 
well as effecting aerosol transport, polydispersisity effects and gas phase species concentration. 

A schematic of aerosol dynamic mechanisms is shown here. The process can begin either in 
the gas phase, through seed particles, or through generation via chemical reactions. Gas phase 
species can also nucleate to form aerosols. These aerosols can then grow through condensation 
of gas phase species on existing particles or trough coagulation of existing particles. 
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Several simplifying assumptions were made to make the model more manageable without 
losing a great deal of generality. 
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k*/=number of molecules in a critical cluster 
ns = number of molecules in cluster at saturation 
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M k /p = e Wt -l or w k = In (M, /p + 1) 




VT = V/ J + V k ' - 



In several of the cases gradients in Mk were very large. These gradients sometimes caused 
difficulties in obtaining the the numerical solution of the governing equations in terms of the 
moments Mk due to, essentially, shocks or discontinuities. 

To alleviate this problem the aerosol governing equations were recast in terms of the a 
transformed variable Wk which is related to the original variable Mk through the mapping Mk/r- 
exp(Wk)-l. This transformed formulation has several advantages over the original formulation 
namely, the range of Wk is smaller (typically ranging from 0 to 50 as opposed to 0 to 1x10 ??) 
reducing the magnitude of gradients. As well, Wk is a conserved variable. Thus, if there are no 
particle dynamics or transport (other than convection) present, Wk is constant, even in a 
compressible flow. The main disadvantage is the the increased complexity (and non-linearity) 
or the RHS source terms. 


